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We study theoretically the Josephson effect and pairing correlations in planar SF1F2S junctions 
that consist of conventional superconductors (S) connected by two metallic monodomain ferromag- 
nets (Fi and F2) with transparent interfaces. We obtain both spin-singlet and -triplet pair ampli- 
tudes and the Josephson current-phase relations for arbitrary orientation of the magnetizations using 
the self-consistent solutions of Eilenberger equations in the clean limit and for a moderate disorder 
in ferromagnets. We find that the long-range spin-triplet correlations cannot prevail in symmetric 
junctions with equal ferromagnetic layers. Surprisingly, the long-range spin-triplet correlations give 
the dominant second harmonic in the Josephson current-phase relation of highly asymmetric SF1F2S 
junctions. The effect is robust against moderate disorder and variations in the layers thickness and 
exchange energy of ferromagnets. 

PACS numbers: 74.45.+C, 74.50.+r 



I. INTRODUCTION 

Spin-triplet superconducting correlations induced in 
heterostructures comprised of superconductors with the 
usual spin-singlet pairing and inhomogeneous ferromag- 
nets have attracted considerable attention recently. 1 
Triplet pairing which is odd in frequency was envisaged 
a long ago^ in an attempt to describe the A phase of su- 
perfluid 3 He. Even though it was found that the pairing 
in superfluid 3 He is odd in space (p-wave) rather than in 
time, it is predicted that even in space (s-wave) and odd 
in time (odd in frequency) pairing does occur in certain 
superconductor (S) - ferromagnet (F) structures with 
inhomogeneous magnetization^! 4 - As a result, supercon- 
ducting correlations can have a long-range propagation 
from SF interfaces, with penetration lengths up to 1/im 
and a nonvanishing Josephson supercurrent through very 
strong ferromagnets 

The first evidence of long-range S-F proximity effect 
came from the experiments on long wires made of Ho 
conical ferromagnet^ and a fully-spin-polarized Cr02 
halfmetallic ferromagnet j£ in which Josephson supercur- 
rent was measured. These experimental findings were 
subsequently verified using different substrates and su- 
perconducting contacts^ Very recently, the effect has 
also been observed in single-crystal ferromagnetic Co 
nanowires^ a Heusler alloy Cr^MnAl^ and synthetic 
antiferromagnets with no net magnetization where thin 
PdNi or Ho layers are combined with Co layers! 10 ' 11 

In SFS Josephson junctions with homogeneous 1 ^ or 
spiral 13 magnetization, the projection of the total spin 
of a pair to the direction of magnetization is conserved 
and only spin-singlet f s and triplet f t o correlations with 
zero spin projection occur. These correlations penetrate 
into the ferromagnet over a short distance determined by 
the exchange energy. For inhomogeneous magnetization, 
odd-frequency triplet correlations fti with nonzero (±1) 
total spin projection are present as well. These corre- 
lations, not suppressed by the exchange interaction, are 
long-ranged and have a dramatic impact on transport 
properties and the Josephson effect^ 



Dominant influence of long-range triplet correlations 
on the Josephson current can be realized in SFS junc- 
tions with magnetically active interfaces ) 14 ' 15 narrow 
domain walls between S and thick F interlayers with 
misaligned magnetizations^ - — or superconductors with 
spin orbit interaction.™ In general, fully developed triplet 
proximity effect can be realized only after inserting a 
singlet-to-triplet "converter," a thin (weak) ferromag- 
netic layer sandwiched between a superconductor and a 
thick (strong) ferromagnet, acting as a "filter" which sup- 
presses the short-range correlations J£r— 

The simplest superconductor-ferromagnet heterostruc- 
tures with inhomogeneous magnetization are F1SF2 and 
SF1F2S junctions with monodomain ferromagnetic lay- 
ers having noncollinear in-plane magnetizations. These 
structures have been studied using Bogoliubov-de Gennes 
equation^ 3 - - — and within the quasiclassical approxima- 
tion in diffusive^ 7 . - — and clea n 35 ' 36 limits using Usadel 
and Elienberger equations, respectively. It has been ex- 
pected that the critical supercurrent in SF1F2S junctions 
has a nonmonotonic dependence on angle between mag- 
netizations due to the long-ranged spin-triplet correla- 
tions. However, in symmetric junctions a monotonic de- 
pendence has been found both in the clean and dirty 
limits, except for nonmonotonicity caused by — it tran- 
sitions^ 3 .^ 

Apparently, the long-range Josephson effect is not fea- 
sible in the junctions with only two F layers ; 17 ' 19 except 
in highly asymmetric SF1F2S junctions at low temper- 
atures, as we will show here. In this case, the long- 
range spin-triplet effect manifests itself as a large sec- 
ond harmonic (1% 3> I\) in the spectral decomposi- 
tion of the Josephson current-phase relation, I(<fi) — 

I\ sin(0) 4- J 2 sin(20) H . The ground state in Josephson 

junctions with ferromagnet can be or 7r stated The en- 
ergy of the junction is proportional to JJ I((f> )d(f> , hence 
a second harmonic leads to degenerate ground states at 
= and <p — 7r. Small contribution of the first har- 
monic lifts the degeneracy which results in coexistence of 
stable and metastable and it states^ 7 - 

In this article, we study the Josephson effect and influ- 
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ence of odd- frequency spin-triplet superconducting corre- 
lations in clean and moderately disordered SF1T2S junc- 
tions with transparent interfaces. Magnetic inter layer is 
composed of two monodomain ferromagnets with arbi- 
trary orientation of in-plane magnetizations. We cal- 
culate pair amplitudes and the Josephson current from 
the self-consistent solutions of the Eilenberger equations. 
Pair amplitudes f s and f t o are short-ranged, while fti is 
long-ranged. 

We show that the influence of misalignment of mag- 
netizations on the Josephson current in symmetric SFFS 
junctions with equal ferromagnetic layers cannot be at- 
tributed to the emergence of long-range spin-triplet cor- 
relations: For thin F layers (compared to the ferromag- 
netic exchange length) all pair amplitudes are equally 
large, while for thick F layers the long-range triplet com- 
ponent is very small. This explains why no substantial 
impact of long-range spin-triplet superconducting corre- 
lations on the Josephson current has been found previ- 
ously for symmetric SFFS junctions both in the ballistic 
and diffusive regimes i 23 ! 32 

We find that long-range spin-triplet Josephson cur- 
rent can be realized only in highly asymmetric SF1F2S 
junctions composed of particularly thin (weak) and thick 
(strong) ferromagnetic layers with noncollinear magne- 
tizations at low temperatures. In that case the second 
harmonic in the Josephson current-phase relation is dom- 
inant. 

The dominant second harmonic provides more sensi- 
tive quantum interferometers (SQUIDs) with effectively 
two times smaller flux quantum 3 ^ and gives the half- 
integer Shapiro steps. 38 The coexistence of and tt 
ground-state configurations in SQUIDs is potentially use- 
ful for experimental study of the quantum superposition 
of macroscopically distinct states^ii Particularly asym- 
metric SF1F2S junctions also enable a robust realization 
of so-called 0-j unctions 4^ In general, Josephson junc- 
tions with a non-sinusoidal current-phase relation are 
shown to be promising for realization of "silent" phase 
qubitsi 2 -^ 

The reminder of the article is organized as follows. In 
Sec. II, we present the model and the equations that we 
use to calculate the Josephson current and spin-singlet 
and triplet pair amplitudes. In Sec. Ill, we provide nu- 
merical results for different ferromagnetic layers thickness 
and strength and orientation of magnetizations. The con- 
clusion is given in Sec. IV. 



II. MODEL AND FORMALISM 

We consider a simple model of an SF1F2S hetero- 
junction consisting of two conventional (s-wave, spin- 
singlet pairing) superconductors (S) and two uniform 
monodomain ferromagnetic layers (Fi and F2) of thick- 
ness d\ and c?2, with angle a = oli — ot\ between their 
in-plane magnetizations (see Fig. [1]). Interfaces between 
layers are fully transparent and magnetically inactive. 




FIG. 1: Schematics of an SF1F2S heterojunction. The mag- 
netization vectors lie in the y-z plane at angles Qi and 0L2 
with respect to the z-axis. 



Superconductivity is described in the framework of the 
Eilenberger quasiclassical theoryji^ Ferromagnetism is 
modeled by the Stoner model, using an exchange energy 
shift 2h between the spin subbands. Disorder is charac- 
terized by the electron mean free path I = vpr, where 
t is the average time between scattering on impurities, 
and vf is the Fermi velocity assumed to be the same 
everywhere. 

Both the clean and moderately diffusive ferromagnetic 
layers are considered. In the clean limit, the mean free 
path / is larger than the two characteristic lengths: the 
ferromagnetic exchange length £p = hvp/h, and the su- 
perconducting coherence length £5 = Hvf/ttAo, where 
Ao is the bulk superconducting pair potential. For mod- 
erate disorder £f < I < t;s- 

In this model, the Eilenberger Green functions 
9 aa '{x,v x ,uj n ), gl a ,(x,v x ,uj n ), f aa '(x,v x ,uj n ), and 

/ , (x, v x ,uj n ) depend on the Cooper pair center-of-mass 
coordinate x along the junction axis, angle 9 of the 
quasiclassical trajectories with respect to the x axis, 
projection v x = vpCOsO of the Fermi velocity vector, 
and on the Matsubara frequencies uj n = TrksTftn + 1), 
n = 0, ±1, . . .. Spin indices are a =t, 4- 

The Eilenberger equation in particle-hole (g> spin space 
can be written in the compact form 



hv x d x g 



w„T3 ® 1 - - A + h(g)/2T,g = 0, (1) 



with normalization condition g 2 = 1. We indicate by ■ 
and • T • 2x2 and 4x4 matrices, respectively. The brack- 
ets (...) denote angular averaging over the Fermi surface 
(integration over 6), and [ , ] denotes a commutator. The 
quasiclassical Green functions 
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are related to the corresponding Gor'kov-Nambu Green 
functions G = —i(T integrated over energy, 
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-f 3 <g> 1 / tfekG 
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where * = (■ip-f^i,^,^) 71 and £k = fi 2 k 2 /2m-/i. The 
matrix V is given by 



V = 1 <g> Re 



h(x) • cr 



+ 173 ® Im 



h(x) • 5 



(4) 



where the components a x ,a y ,a z of the vector a, and 
T\ , 7~2 1 T3 are the Pauli matrices in the spin and the 
particle-hole space, respectively. The in-plane (y-z) mag- 
netizations of the neighboring F layers are not collinear 
in general, and form angles a\ and a 2 with respect 
to the z-axis in the left (Fi) and the right (F2) ferro- 
magnets. The exchange field in ferromagnetic layers is 
h(x) = hi(0, sinai, cos a-y) and h 2 (0, sin 02, cos a 2 ). 
We assume the superconductors are identical, with 



A = 



<t 2 A 
0-2 A* 



(5) 



for x < —d\ and x > d 2 . The self-consistency condition 
for the pair potential A = A(x) is given by 



A = -i\2irN(0)k B Tj2(fn), 



(6) 



where A is the coupling constant, N(0) — mkp /2ir 2 h 2 
is the density of states per spin projection at the Fermi 
level Ep = h 2 k F /2m, and kp = mvp/his the Fermi wave 
number. In F layers A = 0. 

The supercurrent is obtained from the normal Green 
function through the following expression 

I(<t>)=ireN(0)Sk B Tj2J2(v^9<T<T), (?) 

where 4> is the macroscopic phase difference across the 
junction, and S is the area of the junction. In exam- 
ples, the current is normalized to the resistance Rn = 
2v 2 h/Se 2 k 2 F . 

Pair amplitudes, singlet f s , and triplet f t a and f t \, 
with and ±1 projections of the total spin of a pair, are 
defined in terms of anomalous Green functions 

f s (x,t) = -tnN(0)k B Tj2(fn- hl)^*, (8) 
f t0 (x,t) = -i7riV(0)fc s r^{/ ti + / it ) e - i; ""*, (9) 
f tl (x,t) = -tnN(0)k B TY,(fn + fll)^"*- (10) 

In the following we will characterize singlet amplitudes 
by the zero-time f s = f s (x,0). However, since zero-time 
triplet amplitudes identically vanish in agreement with 
the Pauli principle (with or without self-consistency), 25 
we characterize triplets by auxiliary functions using sum- 
mation over negative frequencies only, 

/< = -nrJV(0)fc fl T^(/ n + / ;t ), (11) 

w„<0 



Note that in previous definitions of triplet pair am- 
plitudes the total spin of a pair is projected on the z- 
axis. Physically it is more reasonable to take directions 
of magnetizations in F layers as the spin quantization 
axes. New triplet amplitudes can be introduced in the 
Fi layer (— d\ < x < 0) by simple rotation 

-ifto = cos(ai)(-i/to) - sin(ai)/ti, (13) 
fa = sin(a 1 )(-ifto)+cos(a 1 )fa- (14) 

The same expressions with a,\ — > a 2 hold for F2 layer, 
< x < d 2 . Auxiliary functions / 4 q and are related 
to / t Q and / t "J in the same way. 

For = 0, amplitudes f s and / t "J are real, while / t g 
is imaginary, and for (f> — tt the opposite is true. For 
< < 7r, all amplitudes are complex. In the examples, 
we normalize the amplitudes to the value of f s in the 
bulk superconductors, 



f sb = 2nN(0)k B TY / 



Ac 



(15) 



Here, in the summation over uj n , the high-frequency 
cutoff of 10Ao is used. The temperature dependence 
of the bulk pair potential A is given by A (T) = 

A (0) tanh [\.l^T c jT - l) 4^ 

We consider only transparent interfaces and use con- 
tinuity of the Green functions as the boundary condi- 
tions. For planar SFj^S junctions (three dimensional 
case) Eq. JT]) is solved iteratively together with the self- 
consistency condition, Eq. The averaging over 6 is 
given by (...) = (1/2) JJ(...)d(cos 6). 

Numerical computation is carried out using the col- 
location method. Iterations are performed until self- 
consistency is reached, starting from the stepwise approx- 
imation for the pair potential 



A = An 



e~** /2 G(-x - di) + e^ /2 e(a; - d 2 )] , (16) 



where Q(x) is the Heaviside step function. For finite 
electron mean free path in ferromagnets, the iterative 
procedure starts from the clean limit. We choose the ap- 
propriate boundary conditions in superconductors at the 
distance exceeding 2£s from the SF interfaces. These 
boundary conditions arc given by eliminating the un- 
known constants from the analytical solutions in stepwise 
approximation. To reach self-consistency starting from 
the stepwise A, five to ten iterative steps were sufficient 
for results shown in this article. 



III. RESULTS 

We illustrate our results on SF1F2S planar junctions 
with relatively weak ferromagnets, h/Ep ~ 0.1, and the 



f< = -iTiN(0)k B T 22 (/ft + hi)- ( 12 ) ferromagnetic exchange length £p ~ 20k F . Supercon- 
w „<o ductors are characterized by the bulk pair potential at 
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FIG. 2: Spatial dependence of singlet and triplet pair am- 
plitudes f a , / ( q, and f^, normalized to the bulk singlet 
amplitude We consider a symmetric SF1F2S junction 

(di = d,2 = SOkp 1 ) at low temperature, T/T c = 0.1, with the 
electron mean free path in ferromagnets I = 200fc^ 1 (thick 
curves) and in the clean limit I — > 00 (thin curves). The phase 
difference is <f> = 0. The magnetizations in ferromagnets are 
orthogonal, ot\ — —n/A, 0:2 = f/4, and hi = hi = O.IEf- 
Superconductor is characterized by Ao(0)/i?F = 10 -3 . 



zero temperature A o (0)/Ep = 1CP 3 , which corresponds 
to the superconducting coherence length £s(0) = 636/c^ 1 . 
We assume that all interfaces are fully transparent and 
the Fermi wave numbers in all metals are equal (hp 1 ~ 

lA). 

Detailed analysis is given for symmetric junctions with 
relatively thin [d\ — d,2 = 50k p 1 ) and thick [d\ — d,2 = 
500/s^ 1 ) ferromagnetic layers, Figs.[2H2J an d f° r an highly 
asymmetric junction [d\ = lOkp 1 and c?2 = 990kp 1 ), 
Figs. [5] and [7J In these examples hi = h>2 = O.lEp, 
T/T c = 0.1, and both the clean limit (I — > 00) and mod- 
erate disorder in ferromagnets (I — 200/c^ 1 ) are consid- 
ered. The influences of temperature, thickness and ex- 
change field variations are shown in Figs.|5]and[!|]in the 
clean limit. 

Short-ranged pair amplitudes f s and / t p decay spa- 
tially from the FS interfaces in the same oscillatory man- 
ner. They decay algebraically with length hvp/h in 
the clean limit, and exponentially with the character- 
istic length ^JhD/h in the dirty limit, where diffusion 
coefficient D = vpl/3. 

The long-ranged pair amplitude / t *j is not suppressed 
by the exchange field and penetrates the ferromagnet 
on the scale Hvp/kBT {\Jh~D /ksT) in the clean (dirty) 
limit. In symmetric junctions {d\ = d%, h\ = hi, and 
a,\ = —a/2, a 2 — a /2) all pair amplitudes are practically 
the same in the clean limit (I — > 00) and for moderate dis- 




FIG. 3: The same as in Fig. [2] for ten times thicker ferromag- 
netic films, di = d,2 — 500k p 1 . 
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FIG. 4: Comparison between normalized triplet amplitudes 
2/ti//sf) (solid curve) and rotated 2/ t "j// s i, (dash-dotted 
curve). Parameters are the same as in Fig. [3] for I — > 00. 

order (I = 200fc^ 1 ) in ferromagnets, which is illustrated 
in Figs. [2] and El In symmetric junctions, / t "J does not 
prevail; it can be seen in Fig. [5]that a-dependence of the 
Josephson critical current is monotonic. 

Note that the spatial oscillations of / t "J in Figs. [2]|4]are 
due to our choice of z axis as the spin quantization axis. 
After choosing the magnetization in F layers as the spin 
quantization axis, the rotated amplitude no longer os- 
cillates but instead exhibits monotonic spatial variation 
with the jump at x = (dash-dotted curve in Fig. @|. 
It can be seen that the amplitude / t 5 decays from F1F2 
interface. This is because / t "J in each F layer is gen- 
erated at the interface itself by the projection / t Q sin a 
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FIG. 5: The current-phase relation 1(d)) for a symmetric 
SF1F2S junction with di = d 2 = 500fc^\ hi = h 2 = O.IEf, 
I = 200ftj\ T/T c = 0.1, and for three values of the rela- 
tive angle between magnetizations: a = (solid curve), n/2 
(dashed curve), and n (dash-dotted curve). 



from the neighboring layer. The triplet pair amplitudes 
penetrate into superconductors and monotonically decay 
over the same distance (the bulk superconducting coher- 
ence length £5) as the singlet amplitude saturates to the 
bulk value (see Fig. [4J. 

For thick F layers (when short-ranged amplitude / t Q 
is highly suppressed at the F1F2 interface) the generated 
long-ranged component / t "j is very small. Therefore, in- 
fluence of misalignment of magnetizations on the Joseph- 
son current in symmetric SF1F1S junctions cannot be 
attributed to the emergence of spin-triplet correlations: 
For thin ferromagnetic layers all amplitudes are equally 
large, while for thick ferromagnetic layers the long range 
triplet fti is very small. 

Fully developed long-range spin-triplet proximity effect 
emerges in asymmetric junctions with particularly thin 
and thick F layers (or weak and strong ferromagnets) . 
The thin (weak) F layer acts as a "triplet-generator", 
while the thick (strong) F layer is a "filter" which sup- 
presses the short-ranged components. This is illustrated 
in Fig. [6] for thin Fi layer (d\ — lOfc^ 1 ) and thick F 2 
layer (d 2 = QQOkp 1 ) of equal strength h% = h2 = 0.1 Ep, 
with orthogonal magnetizations, ot\ = —n/2 and a 2 = 0. 
Here, magnetization in the thick (F 2 ) layer is taken along 
the z-axis and along the y-axis in the thin (Fi) layer. 
Because of this choice, the x-dependence of / t "J is mono- 
tonic in F 2 . The thin layer thickness d\ is chosen to give 
maximum triplet current for a moderate disorder in fer- 
romagnets (/ = 200fc i ^ 1 ). This explains why / t "J in Fig. [6] 
is larger for the finite electron mean- free path than in the 
clean limit. 

The current-phase relation for an highly asymmetric 
SFiF 2 S junction (the same parameters as in Fig. [5]) is 
shown in Fig. [7] The critical current for orthogonal 
magnetizations is an order of magnitude larger than for 
collinear magnetizations even though the first harmonic 
is absent. The second harmonic is dominant like at 0- 
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FIG. 6: Spatial dependence of singlet and triplet pair am- 
plitudes f s , / t Q, and f/\, normalized to the bulk singlet 
amplitude f s t- We consider an asymmetric SF1F2S junc- 



tion (di = Wk F and d 2 



at low temperature, 



T/T c = 0.1, with the electron mean free path in ferromag- 
nets I — 200fc^ 1 (thick curves) and in the clean limit I — > 00 
(thin curves). The phase difference is (j> — 0. Magnetization 
in the thin layer is along y-axis, oti = —n/2, and along z- 
axis in the thick layer, a 2 = 0. Here, hi = h% = O.IEf and 
A (0)/£;f = 10" 3 . 
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FIG. 7: The current-phase relation I(d>) for an asymmetric 
SF1F2S junction with di = 10k' 1 and d 2 = 990k' 1 , hi = 
h 2 = 0.1E F , I = 200/c^ 1 , T/T c = 0.1, and for three values 
of the relative angle between magnetizations: a — (solid 
curve), n/2 (dashed curve), and n (dash-dotted curve). 
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FIG. 8: The current-phase relation 7(0) for an asymmetric 
SF1F2S junction in the clean limit I — > 00 with orthogonal 
magnetizations a = ir/2. For di = lOfc^ 1 , d,2 = 990fc^ 1 , 
hi = h,2 = 0.1-Bf three cases are shown: T — 0.1T C (thick 
solid curve); T = 0.01T e (dashed curve); T = 0.9T C (dash- 
dotted curve, magnified 5 x 10 2 times). 




FIG. 9: The current-phase relation I((j>) f° r an asymmetric 
SF1F2S junction in the clean limit I —¥ 00 with orthogonal 
magnetizations a — 7r/2, and T/T c = 0.1. For di = lOfc^ 1 , 
d,2 = 990fc^ 1 , three cases are shown: h\ = hi = O.IEf (thick 
solid curve); hi — O.IEf, hi = 0.2Ef (dashed curve); hi = 
0.05-Ef, h2 = O.IEf (dash-dotted curve). Influence of F layer 
thickness is shown for di — 20k p 1 , d,2 = 980/c^ 1 , and hi = 
hi — O.IEf (thin solid curve). 



7r transitions ) 3 1 except being larger in magnitude and 
long-ranged in the present case. 

Robustness of the long-ranged second harmonic is il- 
lustrated in Figs. [5] and [!|] in the clean limit. The magni- 
tude of the second harmonic is enhanced by lowering the 
temperature, as can be seen from Fig. [8l This is in con- 
trast to the temperature-induced 0-7T transition in which 
the first harmonic is restored away from the transition 
temperature^ 6 - - — When the temperature is close to the 
critical temperature T c , the triplet proximity effect, i.e. 
the second harmonic in 7(0), is no longer dominant in 
an agreement with the results of Ref. Il7l . Small varia- 
tions of the layers thickness and the exchange energy of 



ferromagnets (Fig. [9j, as well as moderate disorder (Fig. 
[7]), do not affect the dominance of the long-ranged second 
harmonic at low temperatures. Since the influence of fcr- 
romagnet is determined by = (h/ Ep)kpd, the shown 
results remain the same for Oi ~ 1 in Fi and 82 ~ 100 
in F2 layers. 

In all calculations we have assumed transparent SF 
interfaces and clean or moderately disordered ferromag- 
nets. Finite transparency much strongly suppresses 
higher harmonics than the first one and the considered 
effect becomes negligible in the tunnel limit. We expect 
the similar influence of large disorder: In diffusive SFS 
junctions at 0-7T transition ; ' 8 where the second har- 
monic is dominant, the Josephson current drops to zero. 
However, this is not the case in the clean SFS junctions 
where the Josephson current is smaller but finite at the 
0-7T transition. 46 



IV. CONCLUSION 

We have studied the Josephson effect in SFiF 2 S planar- 
junctions made of conventional superconductors and two 
monodomain ferromagnetic layers with arbitrary thick- 
ness, strength, and angle between in-plane magnetiza- 
tions. We carried out a detailed analysis of spin-singlet 
and -triplet pairing correlations and the Josephson cur- 
rent in the clean limit and for moderate disorder in fer- 
romagnets by solving self-consistently the Eilenberger 
equations. While the spin-singlet and -triplet correla- 
tions with zero spin projection are short range, the triplet 
correlations with a nonzero spin projection are long range 
and may have a dramatic impact on transport properties 
and the Josephson effect. 

In symmetric SFFS junctions with ferromagnetic layers 
of equal strength, the long-range spin-triplet correlations 
have no substantial impact on the Josephson current both 
in the ballistic and moderately diffusive regimes: For thin 
(weak) ferromagnetic layers all amplitudes are equally 
large, while for thick (strong) layers the long-range triplet 
amplitude is very small. This explains the previous re- 
sults of Refs. m m. 

We have found that fully developed long-range spin- 
triplet proximity effect occurs in highly asymmetric 
SF1F2S junctions at low temperatures and manifests it- 
self as a dominant second harmonic in the Josephson 
current-phase relation. In contrast to the temperature- 
induced 0-7T transition in which the first harmonic is re- 
stored away from the transition temperature, the mag- 
nitude of the second harmonic due to long-range spin- 
triplet correlations increases as the temperature is low- 
ered. The triplet-induced second harmonic is robust 
against moderate disorder and variations in the layers 
thickness and exchange energy of ferromagnets. 

Dominant second harmonic and the resulting ground 
state degeneracy of the Josephson junction (like at 0-7T 
transitions) is experimentally accessible in asymmetric 
SF1F2S junctions with small interface roughness at low 
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temperatures. In addition, the half-periodicity of I(<j>) 
in the considered junctions can be used for quantum in- 
terferometers (SQUIDs) which operate with two times 
smaller flux quantum and can exhibit the superposition 
of macroscopically distinct quantum states even in the 
absence of an external magnetic field.— This has a po- 
tential application in the field of quantum computing. 43 
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